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Systems with long-range interactions display a short-time relaxation towards Quasi Stationary 
States (QSS) whose lifetime increases with the system size. In the paradigmatic Hamiltonian 
Mean-field Model (HMF) out-of-equilibrium phase transitions are predicted and numerically de¬ 
tected which separate homogeneous (zero magnetization) and inhomogeneous (nonzero magnetiza¬ 
tion) QSS. In the former regime, the velocity distribution presents (at least) two large, symmetric, 
bumps, which cannot be self-consistently explained by resorting to the conventional Lynden-Bell 
maximum entropy approach. We propose a generalized maximum entropy scheme which accounts 
for the pseudo-conservation of additional charges, the even momenta of the single particle distri¬ 
bution. These latter are set to the asymptotic values, as estimated by direct integration of the 
underlying Vlasov equation, which formally holds in the thermodynamic limit. Methodologically, 
we operate in the framework of a generalized Gibbs ensemble, as sometimes defined in statistical 
quantum mechanics, which contains an infinite number of conserved charges. The agreement be¬ 
tween theory and simulations is satisfying, both above and below the out of equilibrium transition 
threshold. A precedently unaccessible feature of the QSS, the multiple bumps in the velocity profile, 
is resolved by our new approach. 

PACS numbers: 05.20.-y Classical statistical mechanics; 05.45.-a Nonlinear dynamics and nonlinear dynam¬ 
ical systems. 


Systems subject to long-range forces are known to dis¬ 
play intriguing feature, which manifest both at equilib¬ 
rium and out of equilibrium mm- As an example start¬ 
ing from out-of-equilibrium initial conditions, long-range 
systems get trapped in long lasting Quasi-Stationary 
States (QSS), whose lifetime diverges with the system 
size. Importantly, when performing the mean-field limit 
before the infinite time limit, the system cannot relax 
towards the deputed Boltzmann-Gibbs equilibrium and 
stays permanently confined in the QSS. In this regime, 
the relevant order parameters take values distinct from 
those attained at equilibrium and the system exhibits a 
large gallery of peculiar anomalies, as e.g. non Gaussian 
velocity distributions. QSS have been reported to occur 
for a wide plethora of physical systems for which long- 
range couplings are at play. These include plasma wave 
instabilities [3] relevant for fusion devices, self-gravitating 
systems [1] invoked in the study of non baryonic large 
scale structures formation in the Universe, and Free Elec¬ 
tron Laser mm, light sources of paramount importance 
for their intrinsic flexibility. To elucidate the mechanisms 
which drive the emergence of QSS proved a challenging 
task, that stimulated a vigorous, still open, debate. 

QSS have been interpreted as Vlasov stable station¬ 
ary states, which evolve because of the collisional, fi¬ 
nite N, effects. One can hence gain analytical insight 
into the QSS, by resorting to a Vlasov based maximum 
entropy scheme, as originally pioneered by Lynden-Bell 
in his seminal paper [7]. The aforementioned method 


has been successfully applied both to paradigmatic toy 
models [5] and real world systems [S], allowing to un¬ 
ravel a rich zoology of out-of-equilibrium attributes that 
encompass phase transitions [5] and phase re-entrances 
m- Remarkably, not a single adjustable parameter is 
employed in the theory, which has therefore a fully pre¬ 
dictive value. Despite the general agreement with the 
results of direct simulations, deviations are punctually 
detected when the ergodic hypothesis, which ultimately 
underlies the maximization procedure, fails. Dynamical 
effects need to be properly included into the theory so to 
eventually improve its inherent prognostic ability. Work¬ 
ing along these lines, the so called core halo solution was 
proposed in m- The mechanisms of core-halo forma¬ 
tion shares similarities to the process of evaporative cool¬ 
ing: macroscopic density waves propagates in the hosting 
medium. Particles matching the resonant conditions can 
gain energy at the detriment of collective modes and thus 
leave the inner core to populate a diffusing halo. On the 
other hand, the loss of energy drives a condensation into 
a low energy state of the particles constituting the bulk. 
Due to the Vlasov constraints, the core cannot freeze 
by eventually collapsing into the minimum of the poten¬ 
tial. It in turn approaches the maximum allowed phase 
space density, as imposed by the assigned initial condi¬ 
tions. The method results in a semi-analytical strategy 
to accurately characterized the QSS regime in the regime 
where particles bunch inside the separatrix of a collec¬ 
tive resonance (magnetized phase, as described below). 
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Notwithstanding the undoubtedly success of the theories 
mentioned above, we are however far from having said the 
final word on the genesis of QSS. With reference to the 
celebrated Hamiltonian Mean Field (HMF) model [H], 
the classically recognized arena for challenging QSS stud¬ 
ies, pronounced bumps have been for instance reported 
to spontaneously set in for the homogeneous QSS regime 
mM- These vortices, or phase space clumps, imply 
that the relaxation is incomplete. No explanation of this 
macroscopic, collective phenomenon has been so far pro¬ 
posed. Starting from these premises, the aim of this Let¬ 
ter is to discuss a plausible modification of the original 
Lynden-Bell approach, so to reconcile theory and empir¬ 
ical observations. Although we will develop the analysis 
for the HMF model, the method is general and could 
be in principle mutuated to those settings, from plasma 
physics to Free Electron Laser, where the QSS have been 
shown to occur. 

Let us start by introducing the HMF model, which 
describes the coupled dynamics of N classical rotators. 
The Hamiltonian reads: 

N N 

where Oj stands for the orientation of the j-th rotor and 
Pj is its conjugate momentum. To monitor the evolution 
of the system, one can rely on the magnetization, a global 
order parameter defined as M = |M| = |X;m,|/Af, 
where m.i = (cossinrepresents the local magne¬ 
tization vector. As already recalled, starting from out- 
of-equilibrium initial conditions, the system gets frozen 
in a QSS, whose time duration diverges with the number 
of simulated particles N. When the limit A^ —>■ oo is per¬ 
formed before the infinite time limit (t —> oo) the system 
is stuck in the QSS and cannot proceed towards its associ¬ 
ated equilibrium configuration |15j . In other words, QSSs 
can be interpreted as stationary stable equilibria of the 
system in its mean-field continuum representation. This 
observation immediately translates into a rigorous route 
to analytically inspect the peculiar out-of-equilibrium dy¬ 
namics of systems subject to long-range couplings. 

To this end, we preliminary recall that, in the N ^ oo 
limit, the A^-particle dynamics of Hamiltonian 0 (sim¬ 
ilarly for Hamiltonian of the same class) yields the fol¬ 
lowing Vlasov equation 


dt ^89 dO dp 


( 2 ) 


where f(9,p,t) denotes the microscopic one-particle 
distribution function and 

V{9){f] = 1 - M4f] cos{0) - My[f] sin(d) , (3) 

/ TT pOO 

/ f (9, p,t) cos 0d9dp , (4) 

-TT J —OO 


^4[/] = / / f[9,p,t) s\TL0d9dp . (5) 

J —TT J OO 

The specific energy h[f] = f f(p^/2)f(9,p,t)d9dp — 
(Mf + My — l)/2 and momentum P[f] = 
J J pf{9,p,t)d9dp functionals are conserved quanti¬ 
ties, together with mass normalization. 

The main idea of the Lynden-Bell maximum entropy 
theory, is to coarse-grain the microscopic one-particle dis¬ 
tribution function f{9,p,t). It is then possible to asso¬ 
ciate an entropy to /, the coarse-grained version of the 
single-particle distribution. The sought statistical equi¬ 
librium can be determined by maximizing such an en¬ 
tropy while imposing the conservation of the Vlasov dy¬ 
namical invariants [3 Haim- 

We shall here assume that the initial distribution is 
of the water-bag type. We will in particular restrict the 
analysis to rectangular water-bag in the {9,p) plane, cen¬ 
tered in the origin and of respective half widths Ag and 
Ap. The distribution / takes therefore just two values: 
fo = l/(4AeAp), inside the boundaries of the rectangle, 
and zero outside. The water-bag geometry is entirely 
specified once the energy h[f] = e, and the initial mag¬ 
netization Mq = {MxQ,MyQ) are assigned. 

The Vlasov time evolution can reshape the boundaries 
of the water-bag, while preserving the area inside it. The 
distribution stays therefore two-levels (0, fo) as time pro¬ 
gresses. In this two-level representation, the mixing en¬ 
tropy per particle associated to / reads 

s{f) = - Jdpd9 

(6) 

as it follows from a straightforward combinatorial 
analysis [71 1 16j. The Lynden-Bell recipe amounts to solve 
the following constrained variational problem 


^ln^-h(l-^)ln(l-^ 


fo fo 


fo 


fo 


S{e, cr) = max (s(/) 


Hf) = e-Pif) = a- d9dpf = 1 .(7) 


To carry out the calculations one needs to introduce three 
Lagrange multipliers P/ fo, A//o and fJ./fo for respectively 
energy, momentum and normalization. This leads to the 
following analytical form of the distribution 

g-/3(P^/2-M,,[/] sine-M,,[/] cos e)-\p-fi 
f{^, p) ^ g-/3(p2/2-M„[/] sinO-M^lf] cos ’ 

( 8 ) 

which differs from the Boltzmann-Gibbs expression be¬ 
cause of the “Fermionic” denominator that originates 
from the specific nature of the entropy ([^. By insert¬ 
ing expression Q into energy, momentum and normal¬ 
ization constraints and by making use of the definition of 
magnetization, one can obtain a set of implicit equations 
in the unknowns /3, A, p, and My. Such a system 
can be solved, for any supplied values of the energy e 
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and the initial magnetization Mg, to determine the La¬ 
grange multipliers, together with and My, for which 
the Lynden-Bell entropy is extremal. 

Interestingly, the calculation sketched above enabled 
one to detect an out-of-equilibrium phase transition [9] 
in the reference plane (Mg, e): when assigning the initial 
magnetization and decreasing the energy density e, the 
system experiences a switch from a homogeneous (zero 
magnetization) to a non homogeneous (magnetized) 
QSS. Consequently, the plane of parameters (Mg,e) 
can be partitioned into two adjacent zones, respectively 
associated to an ordered and disordered phase. The 
transition line, which separates the aforementioned 
regimes, is a collection of critical points, as determined 
by the Lynden-Bell theory. Both first and second order 
phase transition are identified, which merge together 
into a tricritical point. The details of the transition 
have been more recently revisited through the core-halo 
approach m, so resulting in a refined description of 
the scrutinized phenomenon. Above threshold, when 
the system is predicted to relax towards a homogeneous 
QSS, two symmetric bumps appear in the velocity 
distribution, as recorded via A'-body and Vlasov based 
simulations. This is the imprint of a collective effect 
which yields the formation of two clusters in the {0,p) 
plane. Neither the Lynden-Bell theory nor the core-halo 
semi-analytic approach could provide a plausible account 
for the spontaneous rise of the bumps, an observation 
which remains to date unexplained. The aim of this 
Letter is to discuss an extension of the Lynden-Bell 
maximum entropy scheme, which accommodates for a 
set of additional Vlasov constraints that prove neces¬ 
sary to benchmark theory and simulations. We will 
here exclusively report numerical simulations of the 
Vlasov equation, the continuum counterpart of the 
original A-body Hamiltonian Q. All simulations were 
performed with the vmf90 program |18j . Provided A 
is sufficiently large, discrete A-body simulations and 
their corresponding Vlasov homologue are virtually 
identical m- Incidentally, we will also show that the 
velocity profiles recorded above threshold return a rather 
intricate gallery of possible structures, ranging from two 
to multiple bumps, an observation which contributes to 
significantly enrich the phenomenology so far reported 
in the literature [23]. 

Consider the 2n-th moment of the dis¬ 
tribution function defined as p 2 n = (p^") = 

p^”/(0,p)(id(ip. For obvious symmetry reasons, 
given the chosen family of initial conditions, only even 
moments need to be considered, the odd ones being 
identical to zero. Starting from out of equilibrium 
initial condition of the water-bag type, with specified 
energy and initial magnetization, one can follow the 
time course of (p^”) for arbitrary choices of the even 
integer n. The typical result of a large campaign of 
simulations is depicted in figure [^: after a violent 


relaxation, as it is customarily termed the sudden 
initial evolution, the moments approach stationary 
plateaux characterized by well defined average values. 
In formulae, limt_>ooP 2 n = e 2 n, where e 2 n are scalars 
determined by the simulations. The Vlasov dynamics 
is sampling its equilibrium solution: the modest oscil¬ 
lations displayed by p 2 n around their mean values e 2 n, 
as shown in figure reflect the local wandering of the 
trajectory inside the target basin of attraction. In the 
discrete realm, finite A effects are indeed responsible for 
the subsequent slow evolution towards the asymptotic 
Boltzmann equilibrium. Working in the continuum 
setting implies silencing endogenous finite A corrections 
and consequently preventing the system from leaving 
the collision-less Vlasov equilibria. These are the QSS 
that we aim at characterizing and which should be 
self-consistently coerced to match the average values 
for the even moments p 2 n, as measured from Vlasov 
simulations. Building on this observation, we elaborate 
on a possible modification of the Lynden-Bell maximum 
entropy solution, which is constrained to reproduce the 
even moments of the distribution, treated as additional 
pseudo-conserved quantities. The moments themselves 
are not Vlasov invariant, but they contribute as dy¬ 
namical constraints to effectively delimit the portion 
of available phase-space. In the following, for purely 
demonstrative purposes, we will make use of the average 
values for (p"), as resulting from the simulations. In 
general it is however possible to construct dedicated 
moments closure scheme which enables to estimate 
the sought quantities, so yielding to a fully predictive 
approach to the QSS characterization [20] ■ 

Mathematically, we shall study the generalized varia¬ 
tional problem 


S{e,a,e 2 n) = max(s(/) 


h{f) = e-P{f) = a- 
dddp/ = 1; P2n = e2n) ■ (9) 


By introducing the Lagrange multipliers i^ 2 n//o to ac¬ 
count for the additional imposed constraints and per¬ 
forming a straightforward calculation one end up with the 
following expression for the extremal distribution func¬ 
tion fQss{0,p) = fi0,p)- 


fgss = fo 


g-/3(p^/2-M,,[/] sine-M,,[/] cosS)-^^ 
^-l3{p^/2-Mylf] sin0-MM] 


( 10 ) 

Notice that we have set A = 0, as it immediately follows 
for symmetry reason. The n-|-2 Lagrange multipliers, to¬ 
gether with the unknown magnetization amount, can be 
determined by imposing the selected constraints. These 
are the conserved quantities, energy and normalization. 
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FIG. 1: Upper panel: time evolution of the first few moments 
P2n with n = 1,2,3. After a violent relaxation the moments 
converges to their asymptotic values £ 2 n. The results refer to a 
choice of the parameter which yields homogeneous QSS. Sim¬ 
ilar observation holds for the magnetized QSS phase. Lower 
panel: running averages in the time interval (250,1000). 


and the stationary values for the even moments of the 
distribution: 



/ fgssdOdp = 1 

^ —OO 

/■°° n - 1 

J P fqssdOdp = e H-^— 

/ OO 

P^'^fqssdOdp = e 2 „ Vn > 2 

-OO 

/ OO 

cos{0) fqssdOdp = 

-OO 

poo 

/ sm{0)fQssd0dp = My 


( 11 ) 


where we recall that fqss depends self-consistently on 
(/3, p, V2n, Mx, My) as dictated by relation (12). In agree¬ 
ment with the original Lynden-Bell scheme, the general¬ 
ized conditions predicts a parameter space (Mo,e) 
partitioned into two region, corresponding to magnetized 
and homogeneous QSS. Above threshold, when the QSS 
magnetization is about zero, equation (121 reduces to the 
following compact expression, only function of p (and so 
called velocity distribution): 


fqssip) /o a 2 /o iri ’ (1^) 

which can be interpreted as an asymptotic free theory. 
The force field acting during the violent relaxation phase 


is implicated in the short time modulation of (p^”), and 
indirectly enter the picture through the quantities e 2 n- 
In this respect, the proposed approach to the QSS re¬ 
sembles a quantum quench from interacting to free the¬ 
ory, as for instance discussed in m, or, more recently, 
in [22j where the inspected system is shown to equili¬ 
brate to a Generalised Gibbs Ensemble [53] . In figure 
the velocity distribution f{p) = f f{0,p)d0 as recorded 
from Vlasov based simulation is plotted (red solid line) 
in the region of a homogeneous QSS, for different values 
of the initial magnetization. The obtained profiles dis¬ 
play the two characteristic bumps to which we alluded 
above. Other collective structures are however present 
and become more evident when Mq is made larger, for 
fixed energy amount (right most panel). 




FIG. 2: The velocity distribution profiles are displayed for 
e = 0.9 and different choices of the initial magnetization (from 
left to right, Mq = 0.3, 0.5, 0.7). In all cases the system is pre¬ 
dicted to converge to a homogeneous QSS. The result of the 
Vlasov simulations are depicted with a (red online) solid line. 
The standard Lynden-Bell theory yields the (green online) 
dashed curve. The modified maximum entropy scheme with 
the inclusion of the additional pseudo conserved quantities 
yields the profiles represented with (blue online) solid lines. 


The velocity distribution predicted by the standard 
Lynden-Bell theory (dashed green) fails to capture the 
fine details of the simulated profile. At variance, the 
generalized theory here described, see eq. (12), which 
sets the even average momenta to the asymptotic values 
determined by Vlasov dynamics, proves definitely more 
adequate in explaining the observations (solid blue). The 
macroscopic details of the velocity profile are adequately 
captured upon truncation at the tenth order (n = 5) 
in the hierarchy of pseudo conserved momenta. No¬ 
tice that the improvements are evident not only in the 
bulk of the distribution, where the resonant particles are 
trapped, but also in the tails, the theoretical curves ad¬ 
hering better to the corresponding numerically computed 
lines. Smaller coherent structures appear in the numer¬ 
ics which are not captured by the theory, at this level of 
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FIG. 3: The velocity distribution profiles are plotted for 
e = 0.6 and different choices of the initial magnetization (from 
left to right, Mq = 0.3, 0.5, 0.7). The system converges to a 
magnetized QSS. The result of the Vlasov simulations are 
drawn with a (red online) solid line. The standard Lynden- 
Bell theory yields to the (green online) dashed curve. Predic¬ 
tions based on modified maximum entropy scheme with the 
inclusion of the additional pseudo conserved quantities are de¬ 
picted with (blue online) solid lines. These latter are almost 
identical to the simulated profiles on the first two panels, from 
left to right. 

resolution. By accounting for higher moments (n > 5) or 
additional pseudo conserved charges might contribute to 
ameliorate the agreement even further. The inclusion of 
additional dynamical constraints is also beneficial in the 
region of magnetized QSS, as it can be appreciated by 
visual inspection of figure ([^ . 

In conclusion we have here discussed a generalization of 
the maximum entropy principle due to Lynden-Bell and 
focus on explaining the emergence of Quasi Stationary 
States in systems subject to long-range forces. By impos¬ 
ing additional dynamical constraints, stemming from the 
Vlasov equation that rules the evolution of the system in 
the continuum limit, allowed us to considerably enhance 
the predictive ability of the theory. In particular, by con¬ 
straining the asymptotic values of the even moments of 
the velocity distribution, we theoretically demonstrated 
that the existence of two (or multiple) islands in the sin¬ 
gle particle phase-space is compatible with a maximum 
entropy principle. Our conclusion is twofold. On the 
one side, we confirmed that a statistical mechanics treat¬ 
ment based on the governing Vlasov equation can be 
successfully invoked to address QSS peculiarities. Dis¬ 
crepancies between empirical observation and the stan¬ 
dard Lynden-Bell theory can be significantly reduced, by 
taking explicit note of key dynamical restrictions. On 
the other side, we expect that our results could open up 
novel avenues to investigate those systems, of both theo¬ 
retical or experimental relevance, where long range forces 
are active and QSSs have been observed (e.g. plasmas. 
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